Comments to the calculation of transverse beam spin asymmetry 
for electron proton elastic scattering 

E. A. Kuraev, S. Bakmaev, V. V. Bytev, and Yu. M. Bystritskiy 
JINR-BLTP, 141980 Dubna, Moscow region, Russian Federation 

E. Tomasi-Gustafsson 
DAPNIA/SPhN, CEA/Saclay, 91191 Gif-sur-Yvette Cedex, France 

(Dated: February 3, 2008) 

Abstract 

The transverse beam spin induced asymmetry is calculated for the scattering of transversally 
polarized electrons on a proton target within a realistic model. Such asymmetry is due to the 
interference between the Born amplitude and the imaginary part of two photon exchange am- 
plitude. In particular, the contribution of non-excited hadron state (elastic) to the two photon 
amplitude is calculated. The elastic contribution requires infrared divergences regularization and 
can be expressed in terms of numerical integrals of the target form factor. The inelastic chan- 
nel corresponding to the one pion hadronic state contribution is enhanced by squared logarithmic 
terms. We show that the ratio of elastic over inelastic channel is of the order of 0.3 and cannot 
be ignored. Enhancement effects due to the decreasing of form factors bring the transverse beam 
asymmetry to values as large as 10 -4 for particular kinematical conditions. 

PACS numbers: 



I. INTRODUCTION 



Polarization observables are a powerful tool for a precise investigation of the nucleon 
structure. Elastic electron proton scattering is the simplest reaction which gives information 
on the dynamical properties of the nucleon, as electromagnetic form factors (EM FFs), if one 
assumes one-photon exchange to be the underlying reaction mechanism. With the advent 
of high energy, highly polarized electron beams and hadron polarimeters in the GeV range, 
it has become possible to measure polarization observables in high precision experiments up 
to large values of the momentum transfer squared. In particular, in elastic ep scattering, 
the polarization of the proton induced by a longitudinally polarized beam allows to access 
the proton EM FFs ratio, as firstly proposed in [TJ [2] . 

The ratio between the longitudinal and the transverse polarization of the proton in the 
scattering plane has been measured up to a value of the momentum transfer squared, 
Q 2 =—q 2 =—t=5.8 GeV 2 , giving surprising results [3J. It has been found that the electric 
and magnetic distribution in the proton are different, and that the electric FF decreases 
faster with Q 2 than previously assumed. However, measurements based on the Rosenbluth 
separation (unpolarized ep elastic cross section) give contradicting results, and confirm that 
the ratio of electromagnetic proton FFs is compatible with unity. 

The extraction of FFs by polarized and unpolarized experiments is based, in both cases, 
on the same formalism, which assumes one photon exchange. To explain the discrepancy 
between these results, it has been suggested that two photon exchange should be taken into 
account [I]. Although such mechanism is suppressed by a factor of a, due to the steep 
decreasing of the electromagnetic FFs with Q 2 , two-photon exchange where the momentum 
transfer is shared between the two photons, can become important with increasing Q 2 . This 
fact was already indicated in the seventies [5], but it was never experimentally observed. 
Recently, a reanalysis of the existing data in deuteron [B] and in proton [7j did not show 
evidence for such mechanism (more precisely the real part of the interference between I7 
and 27 exchange), in the last case at the level of 1%. Although there are a few calculations 
which show that the box diagram contribution is too small to solve the FF ratio discrepancy 
[HI [9], it is interesting to study the implications of the two photon mechanism as one ex- 
pects a detectable contribution of two-photon exchange in electron hadron scattering, with 
increasing Q 2 . The relative role of the two-photon exchange with respect to the main (one 
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photon) contribution is expected to be even larger for d or 3 He, due to the steep decreasing 
of the electromagnetic FFs. The presence of the two photon exchange amplitude should 
be firstly unambiguously experimentally detected. It would induce an additional structure 
function and more complicated expressions for all the observables. An exact calculation of 
two photon exchange in frame of QED for e/x elastic scattering was done in [12], where it 
was found that its contribution to charge asymmetry is very small (at percent level). This 
statement must be valid also for ep scattering, as e/x scattering can be considered to give an 
upper limit for the proton case. 

Other observables, which vanish in the one-photon approximation, can be more sensitive 
to the presence of this mechanism. For example, the transverse beam spin asymmetry 
(TBSA) for the case of scattering of transversally polarized electrons on a proton target 



(the momenta of the particles are shown in brackets and a is the electron polarization) is 



where dVL is the phase space volume of the scattered electron, |a_i_| is the component of 
electron spin, which is transversal to the electron scattering plane, and A(s, t) is the analyzing 
power. TBSA is function of two Mandelstam variables, the total energy s = (pi+p) 2 — M 2 = 
1pp\ and the momentum transfer t = (pi — p[) 2 = —2pip[. Such polarization observable 
was recently measured in high precision experiments devoted to parity violation in elastic ep 
scattering, which have a sensitivity of 1CT 6 i.e., particle per million (ppm). Indeed, nonzero 
values for this asymmetry were found: A = 15.4 ± 5.4 ppm at Q 2 = 0.1 (GeV/c) 2 in the 
SAMPLE experiment, at large scattering angle [UJ whereas in the MAMI experiment A = 



(-8.59±0.89 sia4 ±0.75 s ^) ppm at Q 2 = 0.106 (GeV/c) 2 , and A = (-8.52±2.31 st(lt ±0.87 syst ) 
ppm at Q 2 = 0.230 (GeV/c) 2 [TJ. 



Such observable is sensitive to the imaginary part of the two-photon amplitude, more 
exactly to the interference of Born and imaginary part of box amplitude and vanishes in the 
Born approximation. 

TBSA was recently calculated in in a series of papers. In Ref. [15] a general analysis 
of the double Compton amplitude was derived in terms of 18 kinematical amplitudes. The 



e(a,Pi)+p(p) -> e(pi) +p(p'), 



(1) 



defined as 




da da 



a±\A(s,t) 



(2) 
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kinematics was restricted to the case when the two virtual photons have the same virtuality. 
The results for TBSA at forward angle, in the energy range 3-j- 45 GeV are negative, of the 
order of few ppm, and are not in contradiction with the experimental data. 

In Ref. [13] attention was paid to inelastic hadron states, and the spin asymmetry was 
shown to be enhanced by double logarithmic terms ~ ln 2 (— t/(m 2 e )) for small t kinematics. 
It was stated that the the contribution to TBSA from the proton intermediate state is 
suppressed by several orders of magnitude. The excitation of A(1232) resonance intermediate 
state was also considered. Indeed, A (1232) and any nucleon resonance can not exist in 
the intermediate state as real particles, and a virtual or off-mass shell A resonance has 
no physical meaning. In our opinion, it is more realistic to consider, instead, a definite 
mechanism, with production of a 7r + N intermediate state, which is the simplest inelastic 
channel. 

Inelastic contributions can be expressed in terms of the total photoproduction cross sec- 
tion. This assumption is however an approximation because the exchanged photons are off 
mass shell and t ^ [13] . 

In the intermediate hadronic state, elastic as well as inelastic contributions should be 
taken into account. In lowest order of perturbation theory (PT) the elastic intermediate state 
contribution suffers from infrared divergences, but contributes to TBSA. When considering 
up down proton spin asymmetry, this kind of contribution was not properly considered in 
the previous literature, starting from the basic references [IB] , as well as in more recent 
works. In the present work, the contribution of elastic hadronic state is taken into account 
in terms of a general eikonal phase, and an IR cancellation procedure (IR regularization) 
is applied. After this procedure, it remains a finite contribution to TBSA, which, at our 
knowledge, was not calculated in previous works. 

The purpose of this paper is to point out shortcomings of the previous work, to drive 
the attention to the necessity of considering correctly the elastic intermediate state and to 
compare its relative contribution to the simplest inelastic process of pion production. 

The present calculation differs from the published ones mainly in two respects. We 
show that after applying the infrared regularization procedure, the contribution of elastic 
intermediate states is not suppressed as it was stated in [17] , it produces a finite contribution 
in electron mass zero limit and can be expressed in terms of target FFs. A definite channel 
of pion production in the intermediate state is calculated, in order to estimate the inelastic 
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FIG. 1: Polarized electron- proton elastic scattering: a - Born approximation, b, c - imaginary part 
of two photon exchange amplitude (6 - elastic intermediate state, c - inelastic intermediate state). 

hadronic states. Off-mass shell photons and kinematics with t ^ are considered. For the 
numerical estimation we take a simple ansatz for nucleon FFs, (dipole approximation for 
the Pauli FF and vanishing Dirac FF) which is justified, in principle, only at small or at 
very large values of Q 2 . Nevertheless, it allows to carry on a simple formalism, and does not 
affect qualitatively the results. The present calculation can be generalized to realistic FFs 
in a straightforward way. 

II. FORMALISM 

The amplitude of elastic electron-proton-scattering, Eq. (TXT) 5 taking into account only 
exchanged virtual photons, can be written as [18J: 




where A is a fictitious photon mass and M\ is finite at the limit A — > 0, c is a real constant, 
depending on FFs and it is not calculated here. 

In the language of Feynman diagrams, the phase factor e 1 ^ can be associated with sin- 
gle hadron intermediate hadronic states, whereas Mi corresponds to inelastic intermediate 
hadronic states such as production of neutral pions. 

The specific structure of the amplitude, Eq. [3j suggests a unique way to calculate TBSA, 
eliminating infrared divergences. The asymmetry indeed arises from the interference of the 
Born amplitude with the imaginary part of the two (or more) photon exchange amplitude in 
the lowest order of perturbation theory. The usual mechanism of infrared (IR) singularities 
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cancellation (taking into account real soft photon emission) does not work here, since the 
amplitude of soft photon emission is real. 

The interference of the Born amplitude with the two photon exchange amplitude where 
the intermediate state is a proton, leads to a result which can be expressed in the form 
A(f>(t) + B. Keeping in mind the general form of this scattering amplitude, Eq. |3j in 
order to eliminate the terms containing the phase 4>(t), a regularization procedure, called 
IR regularization, has to be applied. It consists in dropping terms which contain ln(A/m 2 ). 
This is justified by the cancellation of infrared divergent terms when calculating the matrix 
element squared 

The general form for the contribution to the differential cross section for the case of 
transversally polarized incident electron is related to the factor 

m(ppip' 1 a) = me a ^sP a PiPiO' S — mMEE' sin 9 \a±\, (4) 

where m, M, E, E' = E/[l + §(1 - cos 9)] are the mass of electron, the mass of the proton, 
the energies of the initial and final electron, respectively, 9 is the final electron scattering 
angle in the laboratory frame, and a± is the degree of transverse polarization of the incident 
electron. 

The asymmetry is due to the interference of the Born amplitude with the s-channel 
imaginary part of two-photon exchange amplitude. 



A. Elastic contribution 



The elastic contribution to the asymmetry can be written as: 

2 E spm M b M* el 



with 



Si-el V"* I 71 xT 121 I ' 



^|M 6 | 2 = ( 4 ^-Y8(s 2 + u 2 + 2M 2 t)E 2 (t), 

spin 

2^M 6 M* = — F(f) 7T ^ J S e S el ^F{t x )F{t 2 ), (6) 



spin 1 z 

where q 2 2 = — 2e 2 (l — ci^) and e = s/(2ys + M 2 ) in the elastic case. The integration is 
performed in the center of mass system (CMS) of the initial particles, and e is the energy 



of the electron in this reference frame. The proton vertex is described by such model for 
the Pauli and Dirac FFs: = F(t) = (1 + Q 2 /Ql)~ 2 follows a dipole distribution with 

Qq = 0.71 GeV 2 and F 2 {t) = 0. This is a good approximation at very large or very small t 
values. In our case, it is a simple prescription, which allows to factorize the terms containing 
FFs. The traces S e and S e i have the form 

1. 



S el = -Tr{p' + M) 7/i (p + q x + M) lv {p + M) 7A . 



After performing the angular integration (see Appendix [Bj we obtain: 

-stamMEE' sin 9 



A 



el 



el > 



4(s + M 2 )(s 2 + u 2 + 2MH)F(t) 
Q d = -12a g + (4s + 6t)a ls + 2{s + t)a n + [8M 2 - 2{2s + t)]a v - 4M 2 I, 



where the integrals /, a„, an, a ls , a g are given in the Appendix [A] (see Eq. (A4)) 



(7) 
(8) 

(9) 
(10) 



B. Inelastic contribution 



The contribution of the inelastic channel can be written as: 



A 



a ■ t 



inel 



imi- 1 



e'{de'{ 



647T 2- s 2 +u 2 + 2M 2 t J q 2 q 2 S ^el 7r2[ ^ W _ e „ {1 _ c ^ 2 ' 

s ( m^M 2 



dO^dO'Ks 2 - 2Vs + m 2 e'[) 



c w = cospigo, 



m e < ei < 



1 - 



(11) 



2Vs + M 2 V sVs + M 2 J ' 
where the value of g 2 < 3 can be extracted from the total photon-proton cross section 
a tot ~ 0-1 m t ) > a 9 2 /{^ m V)- The trace has the form: 



Sinei = l -Tr(p' + M) 

x (p" + M) 7, 
where, in the inelastic case: 



75 



p'-q + M 
(p> - g ) 2 - M 2 
p-% + M 



p" + % + M 



(p» + g ) 2 - M 2 



75 



X 



(p-q ) 2 -M 2 



p" + g + M 
75 + 75 (p" + g ) 2 -M 2% 



?1,2 



-ee"(l- (3c lj2 ), 



P 2 



m e {e 



ee 



(p + M) 7A , (12) 



(13) 



In the approximation of quadratic logarithms which is valid when e'{ <C M the resulting 
expression for asymmetry can be factorized (in CMS): 

g 2 amM EE's sin 6 



a as 

inel 



64vr 2 (s 2 + u 2 + 2M 2 t)(s + M 2 )F(t) 



RP, 



(14) 



where R is the contribution of the nucleon-pion block: 

-2M 2 + s + 2t r (2M 2 \ / o \ 

R = V s + I ~ ~ 1 ) 7 + + ^ + 



The coefficients 77, 5, 7, J/,, a„, /1, j are given in Appendix [Bj see Eq. (Bl), and the 
term P is the contribution of the 27-intermediate electron block of amplitude: 



dO"t'{de'[ 1 



irqfq 2 t 



in — - hi — in — - H 

77lg mj 2 771 g 3 



III. RESULTS 



Numerical results are presented in Tables |T] and |TT] for different values of the outgoing 
electron scattering angle and of the energy of the initial electron. 

The ratio A e i/A^ el is given in Table [T] (see Eqs. (9fll4). The elastic contribution is not 



negligible, representing more than 30 % of the inelastic one, for most of the kinematical 
region. At small energy and angle it even exceeds the inelastic contribution. The reason 
is that, at small energies, inelastic channels are suppressed in vicinity of the pion pair 
production threshold. 

The total asymmetry A = A e i + A^ el (in units ppm) is given in [TTJ It increases when the 
energy increases, around 9 = 90° where it can be as large as 10 4 . 

The suppression in the vicinity of 9 = or tt, is driven by the sine term. 

In order to compare the calculation with experimental data, different kinematical con- 
ditions should be considered, at low Q 2 . The experimental values or Ref. [12] are well 
reproduced by the present calculation, taking the following values for the cosine of electron 
scattering angle: a = 0.64 and for the pion coupling constant: g 2 = 1.5. 



IV. CONCLUSION 



We calculated the contribution to the electron beam asymmetry in elastic ep scattering, 
which arises from the two photon exchange amplitude in case of non excited hadron inter- 
mediate state. The present calculation shows that that the elastic contribution can not be 
neglected, and that it can be expressed in terms of hadron FFs. 
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The present results contradict previous statements of a specific suppression of this kind of 
contribution [T31 HH HS1 120], the reason being that this source was not properly considered. 
The elastic intermediate state amplitude suffers from infrared divergences (see Eqs. (1), 
(3)). This fact was overlooked in the previous literature. 

Inelastic hadronic intermediate states are enhanced in comparison with the elastic state 
considered above, and can be expressed in terms of double virtual Compton scattering 
(DVCS) amplitude [15J contrary to statement (20], where it was expressed in terms of pho- 
toproduction cross section. 

Double logarithmic enhancement (DL) of inelastic channel takes place at higher orders of 
the QED coupling constant, what can lead to Sudakov form-factor suppresion type. Rigorous 
investigation of higher orders of QED contributions will be subject of further investigations. 

It must be noted that the results are sensitive to the behavior of nucleon FFs, but do not 
chage the qualitatively results of the paper. The dipole like decreasing with the momentum 
transfer squared plays a crucial role and results in a general enhancement of elastic and 
inelastic channels contributions. 

The choice of FFs taken in the present work, allows a simplified calculation. This ap- 
proach is realistic at small Q 2 and at Q 2 ^> 1 GeV 2 due to the suppression of the Fi form 
factor. At small Q 2 , the accuracy of our calculation, due to this approximation is of the 
order of r. 

Another approximation, which affects the precision of the present results, is related to 
the fact that single-logarithmic terms of the form ln(— t/m 2 ) as well as terms of the order 
(m^/M) 2 were omitted. 

Overall, the precision of the present results due to these approximations can be estimated 
at the level of 10 %. 
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E, GeV 


e = 30° 


e = 60° 


6 = 90° 


e = 120° 


9 = 150° 


1.0 


-1.377 


0.138 


0.292 


0.331 


0.340 


1.5 


-0.252 


0.323 


0.379 


0.392 


0.392 


2.0 


0.123 


0.377 


0.389 


0.394 


0.401 


2.5 


0.261 


0.375 


0.378 


0.381 


0.377 


3.0 


0.327 


0.358 


0.363 


0.359 


0.354 


3.5 


0.349 


0.361 


0.339 


0.338 


0.338 


4.0 


0.348 


0.335 


0.327 


0.314 


0.316 



TABLE I: Ratio of elastic over inelastic terms A e i/A\ 



E, GeV 


9 = 30° 


e = 60° 


9 = 90° 


e = 120° 


= 150° 


1.00 


0.76 


-11.70 


-28.27 


-34.39 


-23.037 


1.50 


-2.03 


-22.13 


-43.24 


-43.96 


-26.26 


2.00 


-4.47 


-35.28 


-58.57 


-53.37 


-29.87 


2.50 


-7.42 


-50.05 


-73.77 


-62.71 


-33.39 


3.00 


-11.119 


-66.95 


-88.68 


-71.13 


-37.69 


3.50 


-15.73 


-85.08 


-104.42 


-79.98 


-41.58 


4.00 


-21.58 


-103.16 


-119.07 


-88.17 


-45.32 



TABLE II: Total asymmetry A = A d + A^ el in units of ppm. 

APPENDIX A: ELASTIC CONTRIBUTION (METHOD OF INTEGRATION) 

The biaxial reference frame is used to parametrize the phase volume dO" with fixed 
directions of initial (pi) and scattered (pi') electrons. 

The angular phase space volume of the intermediate real electron is parametrized as 

i c+ 
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D = 1 — a 2 — c\ — c 2 , + 2ciC 2 a = (c 2 — c_)(c + — c 2 ) > 0, 



where 



c± = aci ± y (1 - cf)(l - a 2 ), 

and Ci = cos#i, c 2 = cos #2, a = cos^c, and Q\ = Pip[", 9 2 = PiPi", and C = p\P\' (in CMS, 
where pi+p — 0). Using Euler substitution (t 2 = (c 2 — c_)/ (c + — c 2 )) the relevant integrals 
are 



c+ cfc 2 



7r: 



c + 



dc-i 



71 



VD(1 - f3c 2 ) [(a - p Ci y + (1 - /3 2 )(1 - c 2 )]> 



c+ c 2 cfc 2 



r d 2 c 
J Vd(i - p lCl )(i 



71 



/3 2 c 2 ) Vd 



ln(^ + lnr^ + 21n^^ + ^ 



i-& 2 



4 



d=(l-p 1 p 2 a) 2 -(l-p 2 1 )(l-p 2 2 ), 



where /3i and (3 2 are real number and < (3\ i2 < 1 . 

For the case of a proton in the intermediate state one can write the following kinematical 
relations in CMS: 



Pi+P = Pi +P =Pi +P =0; 

-A 2 + q\ = (pi -p'lf = -A 2 - 2p 2 (l - d) = -2e 2 (l - aCi ); 



g 2 = t 



with 



-2e 2 (l-a) 



l-{3i = - 2 y~e 2 ^ 
p 1 



Then the following integrals can be calculated: 

dO" 



h 
h 



<ql - A 2 )(g 2 2 - A 2 ) 
i r do" i 



4(s + M 2 )' 



7i J q{ - A 2 



(Al) 



(A2) 



Including a realistic hadron form factor F(t), with normalization -F(O) = 1, after applying 
the IR regularization procedure, one obtains: 



IR 



dO"F(q 2 )F(q 2 



(<7 2 -A 2 )(g 2 2 -A 2 ) 

{J, a„(pi + pi),,, a g g^ v + Ou^u* + p' Xil p lu ) + ai s (p liJL p' lu + pi„pij} , (A3) 
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where the scalar coefficients introduced in Eq. (10) have the form: 
d 2 c F{ q \) - 1 F(ql) - 1 , 



1 

2? 

+ 



i |c 2 - a\ 



Dtt 1 - Ci 1 - c 2 



F{qp - 1 F(t) - 1 
l-c 2 



an 



1 

-ten 



e 2 t 2 



ei 2 cF(g 2 )(F(g 2 ) 



1 — a 
1 

' t^ 2 



F(t) 



1 — a 
1 dc 2 [F(ql) - F(t)} 
i 



1 

t7T 



rf 2 C 



rf 2 cg 2 F(g 2 )(F(g 2 )-i; 

TTV^Dfl-Ci) 



L = In 
• i 



|c 2 - a\ 

„2 



+ F(t) • L 



(A4) 



-t(s + M 2 ) 
dc 2 



i l c 2 - a| 



(q 2 2 F(q 2 2 )-tF(t))+tF(t)-L 



It is easy to see that for any form of FFs (particularly for dipole ) all divergences are canceled 
and the final expressions are free from discontinuities. For our choice of FF we have: 



F^h) ~ 1 -4e 2 [Q 2 + e 2 (l - c 1>2 )] F(t) - 1 -4e 2 [Q 2 + 6 2 (1 - a)} 



(A5) 



l-c 1)2 [Q 2 + 2e 2 (l-c li2 )] 2 ' 1-a [Q 2 + 2e 2 (l - a)] 2 ' 

Let us recall that the following relation holds between a, the cosine of scattering electron 
angle in CMS, and the initial electron energy, e: 1 — a = — t/(2e 2 ) , and that g 2 2 = — 2e 2 (l — 
ci, 2 )- 



APPENDIX B: INELASTIC CONTRIBUTION (USEFUL INTEGRALS) 

In the inelastic case, after calculating the traces, the Schouten identity is used to express 
the loop momenta in the numerator via the denominators (we neglect the terms m 2 /M 2 w 
0.02 compared to those of order unity): 

c = (p' - g ) 2 — M 2 = -2p'q , 
b = (p - go) 2 ~M 2 = -2pq , 
d = (p" + go) 2 - M 2 = 2p"q . 
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Finally, one can write the inelastic contribution in terms of the following integrals, Eq. (15): 

s + M 2 



1 

! — 

b 
b 

c 

h l 

c 



doo- 
duo- 

duo- 
dujb 



duJ \ =h = -^ L = ln M 2 
. c 1 — a 

duo- = — —L + a = It, 
b 2(5 



l-af-\(l-~f) (l-a 2 ) 



s + 2M 2 ' 



-2e 2 f3 

IP L 

due = -2e 2 (5 = j 

J du C -q% = {jp + 6p'Y, 7 = 



L + 2a af 

2 



-16M 4 r(l +t] 
1 



[{2M 2 - t)i - 2M 2 j] , 
[(2M 2 -t)j -2M 2 i] , 



du-qo = {a v p' + f3 v QY, (3 V 



- [2Q 2 - s(M 2 + Q 2 )h] , 

[M 2 + Q 2 + 2M 2 sh] , 

jdw<£ = -nQ», "1 = ^ Q = P + Pi, Q 2 = M 2 + s, 
where the following notation holds for the phase volume element: 



d?q Q d?p" 4 „ 

~2^r 5 {Q ~ P - qo) 



sir dO n 
-duo, duo = . 

47T 



2{s + M 2 ) ' 



(Bl) 



(B2) 
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